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, Structures of p- and sd-shell nuclei are studied with the deformed-basis antisymmetrized 

■ molecular dynamics method using the Gogny DIS and Skyrme SLy7 forces as effective inter- 

actions. By the energy variation with a constraint, energy curves as functions of quadrupole 
deformation parameter /3 are obtained. The energy curves for sd-shell nuclei show struc- 
ture change as a function of j3, and suggest shape coexistence. Nuclear structures in the 
deformed region are discussed focusing on deformations and clustering. It is found that the 
deformations often involve cluster structures. Effective-interaction dependence is also dis- 
cussed comparing the results obtained with the two effective interactions. Although the two 
forces give similar results when cluster structures do not develop clearly, they give different 
energies in the largely deformed region when an a cluster develops. 
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§1. Introduction 

> 

' It is well known that a variety of structures appears in excited states as well 

■ as ground states. Shape coexistence is one of the typical examples that a variety 

of structures coexists in a nucleus. Various deformed structures have been observed 
in low-excited states of nuclei in medium- and heavy-mass regions. Recently, su- 
perdeformed rotational band have been observed in ^^Ar^^ and ^'^Ca.^^ Motivating 
OO ! the observation, microscopic studies for deformed structures in sd-shell region have 

' been progressed nowadays. For instance, Inakura et al. have studied deformations 

in sd- and p/-shell region systematically with the Hartree-Fock method using four 
parameter sets of the Skyrme forces, and suggested superdeformations and hyper- 
deformations in these mass region.'^) Moreover, the Hartree-Fock-Bogoliubov (or 
-BCS) -|- generator coordinate method with respect to quadrupole deformation have 
been carried out in sd-shell region. ^^^^^ Thus, nuclear deformation plays an impor- 
tant role also in the sd-shell nuclei, and searching for exotic deformations such as 
superdeformation in this region is one of hot subjects. 

Another important aspect in the light-mass region is cluster structure. It has 
been known that cluster structures often appear in excited states and also in the 
ground states of some nuclei in p-shell and light sd-shell region. Recently, even 
in heavy sd-shell nuclei such as ^^S^^^ and ^'^Ca,^^)'^^) the cluster aspect has been 
found to be important in low-lying deformed states in full microscopic studies. In 
other words, largely deformed states often involve clustering at least in Z = nuclei 
in this mass region. We should point out here that the clustering in the low-lying 
deformed states is not always well marked with an enhanced intercluster distance. 
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but it has spatial correlation of nucleons forming cluster cores. In this sense, we call 
both the formation of cluster cores and the spatially developed clustering "clustering 
correlations" . 

As mentioned above, both the deformations and the cluster correlations play 
important roles in the light-mass region. In the theoretical studies, it is important 
to investigate clustering correlations in deformed states systematically. However, 
such studies have not been progressed enough. 

Our aim is to study structure of deformed states and roles of clustering there. 
For this aim, we study deformations and clustering correlations in p- and sd-shell 
N = Z = even nuclei and a neutron-rich nucleus ^"^O systematically, and analyze 
structure changes as functions of the quadrupole deformation. We use the framework 
of the deformed-basis AMD,^^) which is suitable to study deformations and cluster 
correlations. The deformed-basis AMD wave function is constructed by a Slater 
determinant of triaxially deformed Gauss' wave packets. In this approach, mean- 
field deformations are described by deformations of wave packets while clustering are 
expressed by localizing of centers of the wave packets. We note that the deformed- 
basis AMD is powerful framework to study deformations and clustering in a unified 
manner as already proved in the previous studies of sd-shell nuclei. ^^•'"^'^^ We perform 
energy variation imposing constraint on quadrupole deformation parameter /3 and 
analyze the energy curves and structure changes as functions of /3, focusing on cluster 
correlation in the deformed states. We use two effective interactions, the Gogny DIS 
(DIS)^*^^ and Skyrme SLy7 (SLy7)^^^ forces. The former is a finite-range interaction 
and the latter is a zero-range interaction. Comparing the results obtained by the 
two forces, the interaction dependences are also discussed. 

The paper is organized as follows: Section [2] gives a brief outline of the deformed- 
basis AMD method and properties of the DIS and SLy7 forces. In Sec. [3l the 
calculated results are shown, and features of deformations in p- and sd-shell nuclei 
are described. The interaction dependences are also discussed. In Sec. HI we analyze 
the deformed structures of the p- and sd-shell nuclei focusing on cluster aspects. 
Section [5] is devoted to a conclusion. 



2.1. Wave function and energy variation with constraint in the deformed-basis AMD 

We use the theoretical framework of the deformed-basis AMD method whose 
wave function is a Slater determinant of triaxially deformed Gauss' wave packets, 



§2. Framework 
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\n) = \p) or \n). (2-le) 

Here, the complex parameters Z,j, which represent the centroids of the Gauss' wave 
packets in phase space, take independent values for each single-particle wave function. 
The width parameters Ux , i^y and are real parameters and take independent values 
for each of the x-, y- and z-directions, but are common for all nucleons. The spin 
part \xi) is parametrized by and (3i and the isospin part |rj) is fixed to \p) (proton) 
or |n) (neutron). The quantities {Zj, a^, i^^, f^, i^z} are variational parameters and 
are optimized by energy variation as explained following. 

A trial wave function in the energy variation with constraints is a parity-projected 
wave function, 

where tt is a parity, and Pr is the parity operator. 
The Hamiltonian is, 

H = K + Vn + Vc-Kg, (2-3) 

where K and Kq are the kinetic energy and the energy of the center of mass motion, 
respectively, and Vn is the effective nucleon-nucleon interaction. The Gogny DIS 
and Skyrme SLy7 forces are used in the present work. The Coulomb force Vq is 
approximated by a sum of seven Gauss' functions. 

We perform energy variation and optimize the variational parameters to find a 
state that minimizes the energy of the system E'^ , 

Here, we add the constraint potential V^nst to the expectation value of Hamiltonian 
H in order to obtain local minimum energy states under a constraint. In this study, 
we employ constraint on the quadrupole deformation parameter /?,^°^ 

14„st = t'cnst(/3-/3)'. (2-5) 

When sufficiently large value is chosen for t^cnst; the resultant values of (3 become /3. 
In present study, other quantities such as quadrupole deformation parameter 7 or 
octupole deformation are not constrained. Therefore, those quantities are optimized 
for a given P value by energy variation. 

The energy variation with the AMD wave function is carried out using the fric- 
tional cooling method. The time evolution equations for the complex variational 
parameters Zj,ai and /3j are 

= -,,|^,(. = l,2,..,A), (2.6) 
where Xi is Zj, ccj or /3j, and those for the real parameters v^, a^nd are 

-dr = -^''T^^ (^ = ^,2/,^)- (2-7) 
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Table I. Properties of infinite nuclear matter. 



Force 


Parameter set 


Po,cq [fni"'='] 


ifoo [MeV] 


mo /m 


Gogny 


DIS 


0.160 


209 


0.67 


Skyrme 


SLy7 


0.158 


230 


0.69 




Sill 


0.145 


356 


0.76 



The quantities and are arbitrary positive real numbers. The energy of the 
system decreases as time progresses, and after a sufficient number of time steps, a 
local minimum energy state is obtained under the given constraint. 

2.2. Effective interactions 

We use the Gogny DIS and Skyrme SLy7 forces in the present study, which are 
widely used interactions in microscopic studies of nuclear structure. The Skyrme 
force is zero-range and a density dependent force, and has been widely used to 
self-consistent nuclear structure calculations due to its numerical simplicity. It is 
allocated to make density functionals. The Gogny force consists of a finite-range 
part with Gaussian forms and a zero-range density dependent part. When the Gogny 
force is used in the Hartree-Fock calculations, single-particle wave functions are often 
expanded to harmonic oscillator basis practically. Although fewer studies using the 
Gogny force have been carried out than those using the Skyrme force because of 
more complicated form, it is confirmed that the Gogny force is a credible effective 
interaction by HFB and AMD^^) calculations. 

Let us show the properties of infinite nuclear matter for the DIS and SLy7 forces 
in Table [H In the table, those for the Skyrme SIII forces are also listed for compari- 
son. In the old versions of Skyrme interaction such as the SIII force, the parameters 
are adjusted to fit the data of some nuclei such as doubly closed-shell nuclei. In the 
progress of nuclear matter studies, matter properties are used in adjusting parame- 
ters of modern interactions as well as the ground state properties. For the Skyrme 
SLy7 force, saturation density /3o,eq, incompressibility Koo, equation of state of neu- 
tron matter and so on are adopted as the input data for parameter adjusting. As 
a result, the SLy7 force gives lower incompressibility K^o and effective mass niQ/m 
than the SIII force. 

§3. Deformed structures in p- and sd-shell region 

3.1. N = Z = even p-shell nuclei 

We have performed an energy variation after a projection to a positive-parity 
state imposing the /3 constraint for N = Z = even p-shell nuclei. Figure [1] shows the 
obtained energy curves as functions of matter quadrupole deformation (3 {(3 energy 
curves). 
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Fig. 1. Energy curves as functions of quadrupole deformation parameter 13 for positive-parity states 
of (a) *Be, (b) ^^C, and (c) using the DIS and SLy7 forces. 




3.1.1. ^Be 

The obtained energy curves for 
positive-parity states in '^Be are 
shown in Fig. [^a). In both cases of 
the DIS and SLy7 forces, each energy 
curve has one minimum correspond- 
ing to the ground state at /? ~ 0.7. In 
order to analyze the intrinsic struc- 
ture at the minima, we give the den- 
sity distribution of the intrinsic wave 
function in Fig. [5J As seen in the fig- 
ure, 2a cluster structure develops in 
the energy minimum state. It is con- 
sistent with the well known 2a cluster 
structure in the ground state of ^Be, 
which has been studied by 2a cluster 
models. ■^^•''■^^^ The distance between two a clusters is approximately 2 fm in the case 
of the DIS force. The binding energy of the ground state in the case of the DIS 



2fm 



Fig. 2. Density distributions at the local mini- 
mum on the energy curves (/3 — 0.75) for *Be 
using the DIS force. Symbols "x" indicate 
the centroids of wave packets. 




Symbols "X" indicate centroids of wave pack- quadrupolc deformation parameter (3 

ets. 

Fig. 4. Projection of the energy curves onto the 
/?-7 plane for ^^^O. 



force is approximately 2 MeV larger than that in the case of the SLy7 force. 
3.1.2. 

The energy curves for positive-parity states in ^^C are shown in Fig. [TJb). In the 
deformed region, an oblate shape is favored in both cases of the DIS and SLy7 forces, 
but the shape at the ground states is different between the DIS and SLy7 forces. In 
the case of the DIS force, the oblately deformed state becomes the energy minimum 
at /? = 0.28, while the SLy7 gives the spherical shape of the energy minimum solution. 
Since ^^C is considered to have the oblate shape in the ground band, the DIS force 
gives more reasonable result for the deformation of ^^C than the SLy7 within the 
mean-field approximation. However, even in the case of the SLy7 force, the energy 
curve is quite soft, where the difference of energies at /3 ~ and 0.5 is only 3 
MeV. Therefore, if such treatment as angular momentum projection (AMP) and 
superposition with respect to (5 are incorporated, the oblate shape wave function 
may be significantly contained in the ground state. In other words, those interactions 
give different property in the mean-field framework but they may give similar results 
in a beyond-mean-field framework such as AMP and superposition. It is necessary 
to use beyond-mean-field approaches in order to discuss the detail structures in ^^C 
where the energy curve is soft along the deformation. In large [5 region, 3a cluster 
structures develop. We do not present the energy curves in very large (3 region as 
(3 > 0.5, because it is difficult to obtain convergent solutions in the region. 

3.1.3. 

The energy curves for positive-parity states in ^^O are shown in Fig. He). The 
two forces give quantitatively similar (3 energy curves, in each of which the spherical 
energy minimum for the ground states and a shoulder around /? = 0.7 exist. In 
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the intrinsic wave function around the shoulders, we find a developed q- C cluster 
structure as shown in Fig. [3] for the density distributions. It may correspond to 
the excited band with the developed a-^^C cluster structure which is considered to 
correspond to the K'^ = 0"*" band built on the J'^ = Oj" state. It is interesting 
that the 7 values around the shoulders are different between the results of the DIS 
and SLy7 forces. As shown in Fig. [H triaxially deformed states are obtained in 
f3 > 0.6 region in the case of the DIS force, while axial symmetric shapes are obtained 
in all region in the case of the SLy7 force. It comes from the shapes of ^^C clusters 
composing the a-^'^C cluster structure. As mentioned before, a ^■^C system has oblate 
and spherical shapes at the energy minimal states in the cases of the DIS and SLy7 
forces, respectively. This trend is found also in the excited states of ^^O with the 
a-^^C cluster. In fact, the deformations of ^^C clusters are oblate and spherical in 
the cases of the DIS and SLy7 forces, respectively. As a result, an axial symmetric 
a-^^C cluster structure is obtained because of spherical ^^C cluster in the case of the 
SLy7 force. In the case of the DIS force, a cluster locate on the edge of the oblate 
^^C cluster, and total system forms the triaxial shape. 

Compared to SHF + BCS assuming reflection symmetry using the Skyrme SLy4 
(SLy4) force, the excitation energy around the shoulder at /3 ~ 0.5-0.7 is approxi- 
mately 13 MeV higher, though energy variation after parity projection are performed 
in present study. It is unnatural because the SHF + BCS assuming reflection symme- 
try can not describe an a-^^C cluster state, which is well known structure in excited 
states and expected to be appear in deformed region. The low excitation energy in 
the SHF -|- BCS may be owing to very strong pairing potential. 

3.2. N = Z = even sd-shell nuclei 

We have obtained (3 energy curves by an energy variation after a projection to 
positive-parity states imposing (5 constraint in sd-shell N = Z = even nuclei. We 
have also performed an energy variation for negative-parity states of ^'^Ne which has 
the well known cluster structure. Figures [5] show the energy curves for the 

N = Z = even sd-shell nuclei except for ^"Ca. The detailed study of structures in 
^°Ca using the DIS and SLy7 forces was done in Ref. 14). 
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a) ^^Ne (positive) (b) ^^Ne (negative) 
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(c) ^^Mg (positive) (d) ^^Si (positive) 
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(e) ^^S (positive) (f) ^^Ar (positive) 
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. 5. Energy curves as functions of quadrupole deformation parameter [3 for (a) and (b) '^"Ne, 
(c) 2^Mg, (d) 2* Si, (e) 3^8, and (f) ^^Ar using the DIS and SLy7 force. Figure (b) are the en- 
ergy curves for negative-parity states, and others axe for positive-paxity states. "MF" indicates 
"mean-field-type structure", and "ND" and "SD" indicate "normal-deformed" and "superde- 
formed", respectively (see text). 
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Fig. 6. Density distributions of (a) the ground state (/3 = 0.37), (b) a deformed state in positive- 
parity state {13 — 0.62), (c) a mean-field-type structure in negative-parity states at /3 = 0.52, 



and (d) an a-^^O cluster structure in negative-parity states at /3 
force. Symbols "x" indicate centroids of wave packets. 



: 0.53 for ^"Ne using the DIS 



3.2.1. 20Ne 




J £—L 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
qiiadrupole deformation parameter p 

Fig. 7. Spin-orbit energy curves for negative- 
parity states of ^"Ne. The curves labeled 
"MF" and "a-^^'O" are for mean-field-type 
and an a-^^O cluster structures, respectively. 



The energy curves for positive- 
parity states of ^'^Ne are shown in 
Fig. [5]^a). The two forces give qualita- 
tively similar curves that have prolate 
minima corresponding to the ground 
states at /3 ~ 0.4. Around the min- 
ima, developed cluster structure is 
not seen, but, with the increase of 
the quadrupole deformation parame- 
ter f5, clustering is developing. 
Density distributions of the minimum 
(/? = 0.37) and largely deformed state 
(/3 = 0.62) in the case of the DIS force 
are shown in Fig.[6^a) and (b), respec- 
tively. The binding energy in the case 
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of the DIS force is larger than those 
in the case of the SLy7 force. Compared to HFB using the DIS force, the binding 
energy of the ground state is approximately 3 MeV larger. More detailed discussions 
of positive-parity states in ^''Ne are presented in Ref. 15) and 16). 

The energy curves for negative-parity states in ^'^Ne are shown in Fig. [5jb). In 
the energy variation for a given constraint value of (5, we have obtained two different 
types of structure. One is the deformed mean-field-type structure with no developed 
cluster ("MF") and the other is cluster structure ("a-^'^0"). In the small /? 

region, the energy of MF state is lower than the state, while, in the largely 

deformed region, that of the state becomes lower. The density distributions at 

the crossing point of the energy curves labeled "MF" and in the case of the 

DIS force are shown in Figs. [6l^c) and (d), respectively. Expectation values of the 
spin-orbit term for cluster structure are much smaller in magnitude than those 

for mean-field-type structures as shown in Fig. [71 because the expectation values of 
the spin-orbit term of the ground states in a and ^^O are almost zero. The mean- 
field- type and cluster structures that appear on the curves should correspond 
to the K'^ = 2~ band built on the J'" = 2^ state and the K'^ = 0~ band built on the 
J'" = I2 state, respectively. Comparing the energies of these states between the 
DIS and SLy7 forces, we have found that the energy of the cluster structure 
in the case of the DIS force is approximately 2 MeV lower than that in the case of 
the SLy7 force, whereas the energies of mean-field-type structures are quite similar 
to each other. It suggests that the SLy7 force gives higher excitation energy of the 
cluster state than that of the DIS force. 

3.2.2. 24]vig 

The energy curves for positive- 




parity states in Mg are shown in 
Fig. [5^c). The two forces give qual- 
itatively similar energy curves. The 
mean-field-type structure ("MF") and 
the developed a-^'^Ne cluster struc- 
ture ("a-^^Ne") are obtained in small 
and large f3 region, respectively. The 
curves for the mean-field- type struc- 
ture is quite similar in both cases of 
the DIS and SLy7 forces, while the 
energies of the a-^'^Ne cluster struc- 
tures for the DIS force are approxi- 
mately 2 MeV lower than those for 
the SLy7 force. The mean- field- type and a-^'^Ne cluster structures have minima are 
at /3 ~ 0.4 and 0.8, respectively. Both are prolate deformation. Excitation energies 
of the local minima in the curves of the a-^^Ne cluster structures for the DIS and 
SLy7 forces are 19.9 and 21.5 MeV, respectively. The density distribution of a-^'^Ne 
cluster structure at the local minimum for the DIS force (/3 = 0.84) is shown in 
Fig. El The a cluster locates on the axial symmetric axis of a ^'^Ne cluster. This 



Fig. 8. Density distribution of a-^^Ne cluster 
structure in ^■*Mg (/? = 0.84) using the DIS 
force. 
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quadrupole deformation parameter P 

Fig. 9. Density distribution of the prolate ND Fig. 10. The spin-orbit energy curves for positive- 
state in ^^Si for the DIS force. Symbols "x" parity states in ^*Si. 
indicate centroids of wave packets. 

a-^'^Ne cluster structure with the large prolate deformation is described by four par- 
ticle excitation from sd-shell to p/-shell in the shell-model representation. This is 
consistent with HFB calculation using the DIS force.^) Compared to HFB using the 
DIS force, the binding energy of the ground state is approximately 3 MeV larger. 

3.2.3. 28si 

The energy curves for positive-parity states in ^^Si are shown in Fig. [5]^d). The 
both forces give quantitatively similar energy curves. The curves for prolate and 
oblate shapes as a function of (3 are obtained. The minimum obtained at /3 ~ 0.3 
on the curves for the oblate shape corresponds to the ground state. Since the curve 
for oblate shape is quite soft against quadrupole deformation, the wave function 
in the wide (3 range from spherical to oblate are expected to mix into the ground 
state wave function of ^^Si. The curve for prolate shape has two local minima at 
(3 ~ 0.4 and 0.8. We call those prolate ND and SD local minima, respectively. It 
is considered that the oblate (ground) and prolate band coexist in 28Si,24) which 
is consistent with the present result. The excitation energies of the ND and SD 
minima are estimated to be 9.5 and 16.1 MeV, respectively, of the relative energy 
to the ground state in the case of the DIS force. The prolate ND state in the case 
of the DIS force is slightly softer against quadrupole deformation than that in the 
case of the SLy7 force, whereas the curves are quantitatively similar around the SD 
minima in both cases of the DIS and SLy7 forces. Around the prolate ND minima, 
^^C-a-^^C-like core structures are seen (Fig. [9|), while no cluster structure develops 
around the minima corresponding to the ground state and the SD minima. The 
expectation values of spin-orbit term are shown in Fig. [TOl It is found that the 
energy contribution of the spin-orbit term is remarkably small in the prolate ND 
state. It indicates that ^^C-like core in this state are not the 0^3/2 sub-shell closure 
but they have the spin saturated structure like a 3a structure. In order to discuss the 
details of the coexisting various deformations in 28Si, it is important to incorporate 
the shape mixing in beyond-mean-field approaches. 

3.2.4. 
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Fig. 11. Density distribution of the SD state 
(/3 = 0.67) in using the DIS force. 



The energy curves for positive- 
parity states in ^^S are shown in 
Fig. [5]^e) . The both forces give quan- 
titatively similar energy curves that 
have the prolate ground states at /? ~ 
0.2 and prolate local minima at /? ~ 
0.7 (the SD local minima), respec- 
tively. Excitation energies of the SD 
local minima are estimated to be 12.7 
MeV and 14.6 MeV for the DIS and 
SLy7 forces, respectively. The density 
distribution around the SD local min- 
imum has a neck structure as shown 
in Fig. [TTl It corresponds to ^^O-^^O 
cluster structure of the SD state in ^^S.^^^ Compared to HFB calculation using the 
DIS,^) the binding energies at the ground state and the SD minimum for the DIS is 
4 MeV and 1 MeV larger, respectively. The excitation energies at the SD minimum 
for the SLy7 are a few MeV higher than those with SHF for the SLy4.^^ The HF + 
BCS using the Skyrme SLy6 (SLy6) force gives no SD minimum. ^-^ 

3.2.5. 3^Ar 

The energy curves for positive- 
parity states in ^^Ar are shown in 
Fig. EH^f). Both the forces give quan- 
titatively similar energy curves that 
have minima corresponding to the 
oblate ground state at /3 ~ 0.2 and 
prolate local minima at /5 ~ 0.5 (the 
SD local minima). The 7 values are 
changing as increasing (3 values as 
shown in Fig. [T2l In the small (5 re- 
gion, the system has an oblate shape. 
Then, across the local maxima at /? ~ 
0.35, the system rapidly changes to 
prolate shape. The excitation ener- 
gies of the SD local minima in the 
cases of the DIS and SLy7 forces are 
9.2 and 10.0 MeV, respectively. The 
(3 values of the SD minima are consistent with the calculations SHF^^ and HFB.'') 
The energies at the SD local minimum in the present calculations for the DIS and 
SLy7 forces are almost consistent with that of the the HFB using the DIS and that of 
the SHF using the SLy4 force. The HF + BCS using the SLy6 force gives a spherical 
ground state and no SD minimum.^) Developed cluster structures do not appear in 
all (3 region. 




0.0 0.2 0.4 0.6 0.8 1.0 
quadrupole deformation parameter P 

Fig. 12. Energy curves projected onto /3-7 plane 
of ^®Ar. 



Deformations and clustering in p- and sd-shell nuclei 



13 





3.3. Neutron-rich nuclei O 

The energy curves for positive-parity states in neutron-rich nuclei ^^O are shown 
in Fig. [TSTa). In the case of the DIS force, prolate and triaxial shapes coexist in 
the region from (3 ~ 0.4 to 0.7, whereas the triaxial shape is not obtained in the 
case of the SLy7 force, as shown in Fig. [TSTb). The energy curves for the prolate 
deformations in the cases of the DIS and SLy7 forces are qualitatively similar. Those 
have spherical ground states and local minima at /? ~ 0.7 (the SD local minima). 
The excitation energies of the SD local minima for the DIS and SLy7 forces are 20.0 
and 18.1 MeV, respectively. Neutron density distributions of the axial symmetric 
state and the triaxial structure for the DIS force at /? ~ 0.65 are shown in Fig. [T^a) 
and (b), respectively, for the DIS force. The neutron density distribution of axial 
symmetric structure has a neck, and it is similar to that of the SD local minimum 
in ^^S (Fig. [TT|) . which have the same neutron number = 16 as that of ^^O. 
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§4. Analysis of clustering aspects 

As mentioned in the previous section, it is found that the well-developed cluster 
structure appears in the largely deformed states of p-shell and sd-shell nuclei. In this 
section, we discuss clustering correlations in the deformed states in more detail. For 
this aim, we introduce ideal cluster wave functions given by the Margenau-Brink- 
type wave functions and take the overlap of the deformed-basis AMD wave functions 
with the cluster wave functions. 

We also discuss the difference of the results between the DIS and the Sly7 forces 
from the point of view of the cluster structures. As discussed before, the DIS and 
SLy7 forces give quite similar energy curves for the mean-field-type structures with 
no cluster structure. This is the case in the "MF" states in ^^Mg and ^^Ar, and those 
for the negative-parity states in ^'^Ne (see Sec. 13. 2p . In contrast to those results, the 
difference between the DIS and SLy7 forces is found in the energies of the deformed 
states with the cluster structures. The DIS force gives lower energy for most of the 
cluster-structure states than those of the SLy7 force. When the cluster structures 
develop, the energy difference tends to become larger. Then, it is natural to expect 
that the interaction dependence of the [5 energy curve in the DIS and SLy7 results 
originates in the clustering effects. Comparing the energies of the ideal cluster wave 
functions for the DIS force with those for the SLy7, we will discuss the interaction 
dependence in relation with clustering. 

In order to discuss the clustering effects in the deformed states, we calcu- 
late 2a, equilateral-triangular 3a, a-^^0 and ^^O-^^O cluster structure components 
S2a-, Ssa, •Sq-isq S'i6o-i60) respectively, in the deformed-basis AMD wave func- 
tions for ®Be, ^^C, ^"Ne and ^^S, respectively, as functions of quadrupole deformation 
parameter /3 using Margenau-Brink-type wave functions, which are defined in Ap- 
pendix. The cluster components Sx, where X = 2a, 3a, a-^^0 and ^^O-^^O, take 
< Sx < 1, and, when Sx = 1, it means the deformed-basis wave function is de- 
scribed by a traditional cluster model wave function. We also calculate the energies 
of those cluster wave functions as functions of the intercluster distance using the 
DIS and SLy7 forces. 

Figure [15] shows 2a cluster structure components S2a in the deformed states of 
^Be. In the deformed states of ^Be, the 2a cluster components are largely contained. 
The small S2a value in the small f3 region indicates the breaking of the 2a cluster. 
With the increase of the deformation /3, the value of S2a increases rapidly in /3 ~ 0.3 
region and shows the cluster formation. At the minimum ((3 ~ 0.7) of the energy 
curve, is 0.8 which indicates the 2a cluster structure of the ground state of ^Be. 
Let us show the energies of the 2a cluster wave functions as functions of the distance 
d between the centers of masses of two a's in Fig. [TUl The width parameter of the a 
clusters is fixed to he u = 0.20 fm~^, which is the optimized value for an a nucleus. 
The binding energies of a are 29.7 and 28.1 MeV for the DIS and the SLy7 forces, 
respectively. The local minimum is located at almost the same intercluster distance 
d ~ 3 fm in both cases of the DIS and SLy7 forces, but the energies in the case of 
the DIS force are lower than those in the case of the SLy7 force in all the region 
of the intercluster distance. It means that the 2a cluster structure is energetically 
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quadrupole deformation parameter /3. The 
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Brink wave functions for positive parity states 
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width parameter v is fixed to be v = 0.20 
fm~^, which is optimized for an a cluster with 
the (Os)'* structure. 



favored more in the case of the DIS force than the SLy7 force. Then the difference 
of the (3 energy curve [Fig. [I{a)] between the DIS and the SLy7 forces is interpreted 
as foUows. With the increase of the deformation (3, 2a cluster structure component 
is enhanced. Once the cluster structure develops, the energies of the deformed states 
are different between the DIS and the SLy7 forces reflecting the energy difference 
of the 2a cluster wave function between two forces. As a result, the energy curve is 
steeper in the DIS force, and the deformation with the developed cluster structure 
is more enhanced in the DIS force than in the SLy7 force. 

In the deformed states of ^^C, triangular 3a cluster structure components S^a 
are contained as shown in Fig. [T71 Here the cluster structure component is calcu- 
lated by using the Margenau-Brink-type wave functions with the equilateral-triangle 
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(a) positive (b) negative 




quadrupole deformation parameter p quadrupole deformation parameter P 



Fig. 19. a- O cluster structure components in deformed structures of Ne obtained by energy 
variations imposing constraint of quadrupole deformation parameter /3. In the cases of (a) 
positive and (b) negative parity states are presented. Width parameters v are fixed to = 
Vy = Vz = 0.16 fm~^. 



configurations. In the spherical region, the 3a cluster structure component is almost 
zero, and the component becomes larger as increasing the values of /?. The S-^a is 
approximately 0.4 in the /3 ~ 0.5 region in both cases of the DIS and SLy7 forces. 
This means that the oblately deformed states of ^^C still contain non-cluster struc- 
ture components as well as the 3a cluster components. In Fig. [TSl the energy of 
the 3a cluster wave functions projected to the positive party states as a function of 
the intercluster distance. Similarly to the case of the 2a cluster wave functions, the 
energy for the DIS force is lower than that for the SLy7 force. We consider that 
this energy difference is one of the reasons why the DIS and SLy7 forces give the 
different shapes of the ground state: the former force gives the oblate ground state 
and the latter shows the spherical one. That is to say that the oblate deformation 
is favored in the case of the DIS force because it has the larger energy gain due to 
the 3a clustering. 

Figures [TW a) and (b) show the a-^^0 cluster component 5'q_i6o in the deformed 
states of ^''Ne with positive and negative parity, respectively. The width parameter 
of the a and the ^^O clusters are fixed to = 0.16 fm~^, which is the optimized value 
for ^^Ne in the case of spherical wave packets. As shown in Fig. [5ja), the energy 
minimum state for the ground state is the deformed state with positive parity at 
/? ~ 0.4. In the case of the DIS force, the a-^^0 cluster structure component is 
contained 50 % in the energy minimum state. On the other hand, in the case of 
the SLy7 force the a-^^0 cluster structure component is strongly suppressed in the 
/? < 0.4 region. The SLy7 force shows the rapid change of S'^-ieo around /? ~ 0.5, 
which indicates the drastic transition from the mean-field-type deformed structure 
into the cluster structure. The significant a-^^0 cluster structure component of the 
ground state obtained with the DIS force is consistent with the fact that the ground 
state of ^^Ne have strong a-^^Q correlations.^''''' The larger cluster structure 
component in the case of the DIS force than the SLy7 force is understood by the 
lower energy of the a-^^0 cluster wave function in the DIS force than the SLy7 
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Fig. 20. Energy curves of Margenau-Brink-type cluster wave functions for (a) positive- and 

(b) negative-parity states in the cases of the DIS and SLy7 forces. Width parameters v are 
fixed to Va: = Vy = Vz = 0.16 fm~^. 



[Fig.lMa)]. 

In the deformed states with negative parity of ^'^Ne, two types of structures are 
obtained: the mean- field- type structure and the cluster structure which are 

labeled "MF" and respectively, as mentioned before. As seen in Fig. [TWbl. 

the a-^^O-like deformed states are dominated by the cluster components. In 

contrast, the MF states have no cluster structure component even in the 

largely deformed states. 

Figures [20ja) and (b) show the energies of the cluster wave functions 

projected to the positive- and negative-parity states, respectively, as functions of the 
distance between a and ^^O clusters in the cases of the DIS and SLy7 forces. In 
both the positive- and negative-parity states, the energies for the DIS force are lower 
than those for the SLy7 force as in the cases of 2a and 3a cluster wave functions. 
Reflecting this energy difference of the cluster wave function between the DIS and 
the SLy7 forces, the energies of the deformed states with clustering obtained by the 
/^-constraint are relatively lower in the case of the DIS than the case of the SLy7 
force. In particular, those for the negative-parity states are almost consistent with 
the energies of the negative-parity cluster wave functions. This results in the lower 
excitation energy of the negative-parity states with the cluster correlation which 
constructs the K'^ = 0~ band in the ^'^Ne. 

Figure [21] shows ^^O-^^O cluster structure components 5i6o-i6o in the deformed 
states of ^^S. The width parameter of the ^^O cluster is fixed iov = 0.16 fm~^, which 
is the optimized value for a ^^O cluster. In the small deformation region (/3 < 0.4), 
the ^^O-^^O cluster structure component is almost zero, while it is significantly 
contained in the largely deformed states(/3 > 0.6). The value of 5i6o-i6o of the 
local minimum for the SD state at /3 ~ 0.7 is 0.3, which indicates the ""^O-^^O 
cluster correlation in the SD state. Figure [22] shows the energy curves projected 
to positive party as a function of the distance between two ^^O clusters for the DIS 
and SLy7 forces. The DIS force gives the lower energy of the cluster wave function 
than the SLy7 forces. It is consistent that the energy of the SD local minimum 
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on the /3 energy curve for the DIS force is shghtly lower than that for the SLy7 
force, and the SD state in case of DIS force contains larger ^^O-^^O cluster structure 
component S'leo-ieo- 

The present results show that the cluster structures affect the deformed states 
in sd-shell nuclei as well as p-shell nuclei. Although the DIS and SLy7 forces give 
the similar energy curves for the mean-field-type structures, they give the different 
energies for the ideal cluster wave functions. In the p-shell and sd-shell regions, the 
DIS force gives the lower energies of the cluster states than those of the SLy7 force. 
This difference is reflected in the excitation energies of the largely deformed states 
with chistcring correlation. In other words, the energy of the cluster structures are 
reflected on the /3 energy curves due to the clustering effect in the deformed states. 

§5. Conclusion 

We have studied deformations and clustering correlations in the p-shell and sd- 
shell regions using the Gogny DIS and Skyrme SLy7 forces. The effective-interaction 
dependence has been also discussed. We have obtained the energy curves as functions 
of a quadrupole deformation parameter /3 with a parity-projected energy variation 
imposing a constraint of quadrupole deformation parameter /3 to the deformcd-basis 
AMD wave functions. We have found that the various deformed structures appear 
with the increase of the deformation /? in both the cases of the DIS and SLy7 forces, 
which give qualitatively similar results. 

We have found that various cluster structures appear in p- and .sc?-shcll regions. 
For example, in the large /3 region of = Z < 12 nuclei, the well developed a cluster 
structures appear in both the cases of the DIS and SLy7 forces, and the ^^O-^^O 
clustering is found in the SD local minimum of the /3 energy curves of ^^S. On the 
other hand, no cluster structure is seen in the deformed states in ^®Ar and ^^O. 

In comparison of the results between the DIS and the SLy7 forces, the properties 
of the /3 energy curves are qualitatively similar except for ^^O, and the 7 values are 
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also similar except for O and O. In O, the two forces give different 7 values in 
the large (5 region because of the different shapes of -^^C clusters. In ^^O, the triaxial 
deformation appears in the /3 region between the spherical shape for the ground state 
and the SD local minimum in the case of the DIS force, whereas the SLy7 force gives 
only the axial symmetric shapes. Cluster structures in the deformed states for the 
DIS and the SLy7 forces are also qualitatively similar to each other in most cases. 

For the mean- field- type structures in the sd-shell region, the two forces give 
quantitatively quite similar results. In contrast, we have found the strong effective 
interaction dependence in p-shell and the deformed states with strong cluster corre- 
lations in sd-shell nuclei. For example, the quantitative difference is found in ^^C, for 
which the DIS force gives the oblate shape while the SLy7 force gives the spherical 
shape at the energy minimum in the /3 energy curve. The cluster structures in the 
deformed states of '^Be, ^^C, ^°Ne and are analyzed using 2a, 3a, a-^^O and 
16q_16q Margenau-Brink-type cluster wave functions, respectively. Although the 
DIS and SLy7 forces give the similar energy curves for the mean-field- type struc- 
tures, they give the different energies for the ideal cluster wave functions. In the 
p-shell and sd-shell regions, the DIS force gives the lower energies of the cluster 
states than those of the SLy7 force. This difference is reflected in the excitation 
energies of the largely deformed states with clustering correlation. 

In conclusion, various cluster structures appear in sd-shell region as well as p- 
shell region. In studying nuclear deformations, cluster structures and properties 
of lighter nuclei should be taken into account, because they can affect the nuclear 
deformations through the cluster correlation in the deformed states. 
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Appendix A 

Definition of a cluster structure component 

In present paper, cluster structure components Sx, where X = 2a, 3a, 
or ^^O-^^O, arc evaluated in order to discuss clustering correlation. In this section, 
we define the cluster structure component Sx in analogy of cluster probability of 
AMD wave function in Refs. 12), 27). 

The normalized vr-parity projected deformed-basis AMD wave function l^'^) is 
generally divided into X cluster component \^x) and the residual part 1^^^^) as 

1$-) = c 1$^) + (A-i) 
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here 

(^xl$R,)=0. (A-2) 

We introduce the projection operator Px which projects out the X cluster component 
from the l^'^), 

Px\^^) = c \^\) = cA\xJ{^{Xi)). (A-3) 

i 

where 4>{Xi) is the normahzed internal wave function of the cluster Xj which is 
contained in X cluster structure and x is the wave function of relative motion of the 
clusters. 

The projection operator Px is defined as following. Suppose the wave functions 
1$^'^) that span the functional space of X cluster structure. The wave function 
is a Margenau-Brink-type wave function with intercluster distance dj, which is 

di = dj^in + ^^"^ (0 < i < imax), (A-4) 

Ad=^^^^^^, (A-5) 
^max 

where d^^-^ and (imax are minimum and maximum values of intercluster distance, 
respectively. In present study, dmin and Ad. arc adopted to (imin ~ 0.1 fm and Ad, ~ 1 
fm, respectively. Orthogonalized wave function l^^'^) is obtained by performing 
unitary transformation from which is vr-parity projected state of 

\^^-) = u^,\^f-), (A-6) 

(ia"l*f")=^a/3. (A-7) 

Using the set of orthonormalized wave functions jij'^), the projection operator Px 
to project out the X cluster cluster component is defined as 

a 

Using the projection operator Px, cluster components Sx is defined as 

Sx = |c|2 = (^>^|Px|$") = ^ (A-9) 

a 

where wave functions |$^) and |$) are rotated so as to diagonalize moment of inertia 
matrices as 

< < (^2^. (A-IO) 
Since l^^'^) are orthonormalized, Sx takes < 5x < 1. 
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